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Abstract 

> _ 

r*"** | Spherically symmetric static fluid sources are endowed with rotation 

■ and embedded in Kerr empty space-time up to an including quadratic 
terms in an angular velocity parameter using Darmois junction condi- 

■ tions. Einstein's equation's for the system are developed in terms of linear 
\ ordinary differential equations. The boundary of the rotating source is ex- 
pressed explicitly in terms of sinusoidal functions of the polar angle which 

| differ somewhat according to whether an equation of state exists between 

O" 1 internal density and supporting pressure. 

(3JT), Following the publication by Kerr [1] of the metric which describes ana- 

lytically, the asymptotically flat, vacuum gravitational field outside a rotating 
source in terms of Einstein's field equations, there has been much discussion 
concerning the existence of possible interior solutions which match the exterior 
smoothly. In an important development Hartle [2] uses a second order pertur- 
bation technique to describe the slow rotation of equilibrium configurations of 
cold stars having constant angular velocity. Solutions of Einstein's equations 
are developed in terms of Legendre polynomials but the issue of matching the re- 
sults to Kerr empty space-time is not addressed. In the case of non-equilibrium 
configurations Kegeles [3] has applied the method to Robertson- Walker dust 
sources up to the first order in angular velocity parameter although, the re- 
sults are somewhat restrictive and are not suitable for application to sources 
supported by internal pressure. In a recent work the case of the Wahlquist [4] 
closed form interior was shown not to fit the Kerr exterior by Bradley et al [5] . 
Only for the important case of thin super-massive rotating discs, supported 
by internal pressure have analytic sources for the Kerr metric been found (Pi- 
chon and Lynden-Bell [6]). This has led to an 'embarrassing hiatus' according 
to Bradley et al [5] in the number of potential interior solutions available 
for matching which in turn has contributed to a lack in the development in 
the theory of differentially rotating fluid bodies in general relativity. Yet it is 
important to develop further the relativistic theory of rotation since it has con- 
siderable potential application in astrophysics, for example, in the description 
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of the gravitational collapse of rotating matter, quasars, or potential sources 
for gravitational radiation. 

It is the aim here to focus upon one aspect of the problem of matching rotat- 
ing interior solutions of Einstein's equations to the Kerr metric. In particular it 
is the intention to show how static perfect fluid sources endowed with rotation 
to an accuracy of second order terms in an angular velocity parameter q may 
be matched to the Kerr metric also expressed to an accuracy of quadratic terms 
in angular velocity parameter. In this sense the method employed is similar to 
that of Bradley et al [5]. Moreover, the results given here apply to fluid inte- 
riors having variable interior angular velocity and so to that extent generalise 
the discussion of Hartle [2]. Indeed the problem considered here is essentially 
that of matching the generalised Hartle rotating interiors to Kerr empty space- 
time although the approach here differs significantly. The matching problem is 
addressed by means of the application of Darmois [7] junction conditions which 
are discussed widely in the literature for example, Misner et al [8], Mars & Sen- 
ovilla [9], Stephani [10], Hernandez-Pastora et al [11] and Bonnor and Vickers 
[12]. In the Darmois approach it is necessary that the components of the metric 
tensor, and also the extrinsic curvature for the Kerr exterior and the interior 
source are continuous at the boundary surface. 

In the following Einstein's equations for a perfect fluid source will be ex- 
pressed in the form: 

G a b = -87tT b a , T b a = (p + p)u a u b -S a bP , (1) 

where p, p are the respective rotating source density and supporting internal 
pressure and u a are the components of the velocity four- vector with the property 
that u a u a = 1. 



1 The fluid source and Darmois junction conditions 
with empty space-time 

Consider the well known spherically static perfect fluid source which is repre- 
sented here by the metric: 

da 2 s = e 2X drf - e^dt 2 - £ 2 d9 2 - f sin 2 9d(p 2 (2) 

where A = A (£) and p = fi (£) such that 

% = -A| + ^ + ^±^-^ + l (3) 

and a suffix indicates a derivative. The respective supporting pressure and 
density of the static cases are: 

8np s = ^ + — - ^ (4) 
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and 



8ttp s 



2e~ 2 ^ 
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(5) 



In the following this source will be endowed with rotation up to and includ- 
ing quadratic terms in angular velocity parameter q. It is assumed that the 
resulting physical system may be described by means of the metric 



da 2 = e 2X (1 + Qq 2 ) dr, 2 - e 2 ^ (l + Uq 2 ) df - f (l + Vq 2 ) dB 7 



-£ 2 (1 + Wq 2 ) sin 2 



2X£ 2 gsin 2 6d<pdr] 



(6) 



where U,V, W and Q are each functions of £ and 9 whilst X is a function of £ 
alone. This particular gauge has been chosen since it is the intention to match 
this interior source to empty space-time using the well known Kerr metric 
expressed in Boyer and Lindquist coordinates at a boundary between source 
and empty space-time given by: 
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where / = / (9) and £ 1S constant. All expressions will be accurate up to and 
including quadratic terms in q so that the Kerr metric may be written as: 

da 2 = ^ AA drf ~ 2^439 sin 2 6>#g?t? - ^ndS, 2 - ^ 2 2<i9 2 - ^ 33 sin 2 9d(f) 2 (8) 



where: 



^44 
*43 
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*22 
^33 



1 

2m 

T 



2 m ( q 2 cos 2 9 



1 



i 2m 



1 + 



£ 2 + q 2 cos 2 9 
2mq 2 sin 2 9 



q 2 cos 2 9 



+ q 2 + t 2 



e i 



2m 



(9) 



Darmois junction conditions are employed at the boundary so that both the 
metric tensor g a h and the extrinsic curvature K a ^ must be continuous on the 
bondary surface (0. Using the suffix b to indicate evaluation at £ = £ fe and the 
suffix to indicate evaluation at £ = £ it is straight forward to show that the 
following conditions on A, fi, U, V, W , Q and X and derivatives with respect 
to £ hold at the boundary: 
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Vb 



Co 

1 -K 



2K£, 

K cos 2 -1 , , , /(^-l) 
cos 2 6 

4/{/»J _ 2cos 2 fl 
P P 3 

1 + (1 - K) sin 2 6 



e 



3(if-l)sin 2 fl-2 4/{ M J 

1^/6 = 73 + 



Qb 



Co Co 
(1 - if) cos 2 fl 



(2if 2 - g - 1) cos 2 g /(1-3K)VJ /(if 2 -l) 
= 737^ + TIF + 



X b ~ — 2 



Co^ 2 Co^ 2£ 2 ^ 2 

1 -K 

cT 



- ^ do) 

So 

For future convience it is useful to define the functions H (£, 0) and h (£) such 
that: 

#(C,0) = Q (£,#) + sin 2 9fX 2 e~ 2X 



so: 



3sin 2 fl(l-if) (2if 2 -if-l) /(if 2 -l) 
1 ~ $ K\l ~ + 2^ 2 

/{mJ (1-3JQ 



£o^ 

fc 6 = 3£ (tf-l) (12) 
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2 Einstein's equations for the rotating source 

There are three non-trivial Einstein equations which need to be solved for 
a fluid source in the gauge © • In the first the condition T\ = = T 2 gives rise 
to: 

-u e x, + H e x, + We ,-^ + He ,-^ + cos 9 ^ = (13) 

£ q sinf 

whilst T-f = T| becomes: 

e~ 2fl f .. . ... . ... 
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cosfl ([fr + y -2W fl ) _ H ee + W 0fl £f - F 

2£ 2 sin# 2i 2 e 

h 2 sin 2 
2£ 6 

= (14) 
The final equation is (T3 + p) (Tf + p) — TfT^ = which is explicitly: 

e~ 2 » (Af + ^) - A^ W { - ^ + + % 

2y 5 I cosfl (H e + U e ){X, + ^) 

| (H ee + Efo) (\ + gg) _ . 2 fl , (Vj^g) fe2 > \ 

+ ^ Sm ^5 + ^6 J 

= (15) 

An important differential consequence of these equations, since it is inde- 
pendent of any second derivatives with respect to £ is the relation: 

^^ { u -H e )e\dk-»d 

+e -2/i _ ^ ) ( £2 A? _ ^2 + 2 ^ ^ + ^) 

e~ 2 ^ A^ sin 9 , n , , s 

(e 2 ^-£C/ e -£ 2 ^ + £^) (A{ + m € ) 



cos # 
sin $ 

+ a ( W eee -2We + 3V e + H m - H 6 ) (A ? + liA 

e _2/1 (At + Lit) sin 6* , , „ „ . 

+ n £ 2 ^V - £ ^ + £ 2 + £ 

cos ' 

COS $ 

-£ + tV ei; +2H e } + — (-2 iy e + y e - (A e + 



+ (2Wee-Vee + Hee){X i + ^)+sm 2 9 , 
(16) 
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The expression for the supporting internal pressure is: 

2u 2 (We\t V £ \ £ 2U\ £ W £ V £ U H £ 

q 2 

( w oe + 2V + H ee) 

q 2 cos 9 (2W e -V e + H ) h 2 q 2 sin 2 9 
+ 2£ 2 sin0 + 4£ 6 1 ° 

and the fluid density is: 



i„ 2 (W £ \ £ V £ X £ W £ 5V £ 
87rp = 8np s + e- 2 V | -LI _ _L_i _ _| _ V ££ _ _A + ^ 

U £ 2p £ U U) q 2 , 
, g 2 cos (-2W* + y e + fl*) , /i 2 g 2 sin 2 

+ 2^rnl + 4£ 6 ^ 

The angular velocity may be written as: 
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T 4 3 _ qh £ e x+ » 



3 Description of the boundary 

Application of the boundary conditions ifTUj) and 1)12(1 show that equation <|13|) 
is satisfied identically whilst (|17|) gives: 

8vrp 6 = (20) 

as expected. 

Whenever \ £ + fj, £ ^ 0, so that from the second of ljTU|) 

J = 2^ {^} Q K-K + 1^0 (21) 

then equation l|T5jl may be used to define / (9) since at the boundary it reduces 
to: 

-2 



2 cos sin 9KJ (4£ {/xJ Q if (K - 1) + 7K 2 - 10K - 1 

feo {mJ ^ + ^ - l) (2^o {mJ K (K + 1) - 3K 2 + 4K - 1 
— cos sin 9K 2 {h 2 } b 
= (22) 

Thus in such cases it follows that: 

/ (9) = a + ai sin 2 9 (23) 
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where ao is an arbitrary constant and a\ is found using (|22[). 

However the case {A^ + /^} = is also of importance. From equation 1)16(1 
it follows that {h^} b = and so the value of the angular velocity L{^ b ,rj) at 
the boundary surface must be determined by careful analysis of the limiting 
behaviour of (|19|) at £ = £ 6 . Moreover since this case corresponds to the con- 
dition J = then the second of (fT2l and © imply that the surface density 
for the static solution must satisfy {8irp s } = 0. This will invariably be the 
case in physical situations where an equation of state of the type p = F (p) is 
additionally imposed on equations ((T3|) to (fT3|) to close the system of Einstein's 
equations. For such cases the condition {8np} b = also gives: 



c3 
so 



{U(\ Q 8cos 2 6> 



, f ( 2K Mo | 
+Jee 1- 



e 



so 



1 + K 

so 



(24) 



This equation may be used to define the equation of the boundary / = / (6) 
whenever {U{\ Q and {V^} are known. Clearly whenever {U^} and {^} 
depend purely on terms in sin 2 6 then / (9) will have a solution of the form 1)23)1 
and ao and a\ would be found using (|24|1. 



4 Generating solutions of Einstein's equations 

Consider now the condition { Ag + ^ and also the condition {A^ + /^} = 
provided that {U^} and {^£§} depend purely on terms in sin 2 0. The case 
when {U^} and {V^} depend on higher harmonics is discussed in the next 
section. 

It follows that / (9) has the form 1)23)1 and for such cases solutions of Ein- 
stein's equations satisfying the boundary Darmois boundary conditions may be 
found by writing: 

U(£,0) = Go (0 + d (0 sin 2 9 
V(Z,9) = G 2 (0 + G 3 (0 sin 2 
= G 2 (£) + G 5 (£) sin 2 
#(£,0) = G 6 (0 + G 7 (0 sin 2 (25) 

The equations ((T^|) . 1)14)1 and lfT5|) then give rise to four independent equa- 
tions for Gj (£) as follows: 

2G 7 A e -2G 1 A 5 -^-^ + 2G 75 +3G 5e -G 3e = (26) 
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9 r \ r \ r \ G it G H G \ G k 
2G 7( A ? - G 35 A 5 - ti e A{ - — + — — 

4G 7 e 2 ^ 8G 5 e 2 ^ 4G 3 e 2fl /i 2 e 2 ^ 
~ G 7(f J '£ ~ G 5 ( ^ + G7 e5 + Gs^ 

(27) 



2G 7 ^ A e - — - ^ Ag - — + —2- - ^ J + 

2Gi f - A? + ^ + /i f A e + ^ + , k t e ^ — r- 

v € e 55 ey 4£ 5 (A € + / ^) 

+2G 7e (A ? - ^) + G h (h + j- 
+G 3s (-3A ? - | + ^ + 



16G 5 e 2 ^ 8G 3 e 2 ^ /i 2 e 2 ^ 



(28) 



+ 



r t 

2(G 7 + 3G 5 - G 3 + G 2 - d - G ) e 2 » 



-2Gq^X^ + G2 f X^ + Gq^A^ H — ^ H — ^ 



c 2 



+ Gfyfj,^ + G^/x^ — G6 e5 — G2 ?e 
= (29) 

The first three contain only G\ (£), G 3 (£), G5 (£) and G 7 (£) whilst the 
fourth contains each of the unknown Gj (£). This under-determined set may be 
closed with the aid of an imposed equation of state as discussed in the previous 
section. 

The boundary conditions (fl7)|) and (|T2|) together with (|25|) then give rise to 
the following: 

Go(e ) = V(rl) +2a °H 

1 ai (1 - if) 



,0 

Gi^o) = 2ai{^} 
C2 (Co) = 71 

G 3 (e ) = -i 

Co 

G 5 (Co) = ^ 
So 

G 6 (Co) = 

G 7 (Co) = ^ (30) 
so 
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and: 



G 2 , (Co) 
Gs, (Co) 
G k (Co) 
Ge £ (Co) 



4{/i ? } a 2 



Co 



r 

so 



4{mJ «i , 2 



Co 



+ 



t 3 
SO 



4{//J ai 3(-KT-l) 
Co 



£ 3 

so 



a 



Co 



+ 



2 Co Co 



G 7 e (Co) = — 7T 



Co 5 Co^ 



«1 



2C^ 2 



ai 3 



c ^ eg* 2e^ 2 c^ 2 

3oo{mJ ( 



3ai{/xJ ( 



(31) 



5 Discussion and conclusion 

In the analysis above it has been shown how the problem of determining slowly 
rotating sources embedded in Kerr empty space-time has reduced to the prob- 
lem of determining solutions of the system of ordinary differential equations 
(f2lj|) to ((2~§]) subject to the conditions (|3"U|) and l|3"Tj) , Whilst these equations 
are under-determined they may be supplemented by an equation of state to 
make unique solution possible. Clearly the solution of the equations is straight 
forward enough although analytic solution in closed form for any known A and 
\x seem not to be possible. 

The most appropriate techniques for solution are either numerical or al- 
ternatively a Taylor's series solution approach for Gi (C) about C = Co- Two 
variations on the series method are worthy of note. In the first one may assume 
that both A and fi are known, for example, the Schwarzschild interior solution 
may be specified, and then take the system (|2f))) to subject to the condi- 
tions (123) and (jHU linear set of equations for Gi (C) ■ In the second one may 
assume that neither A and [i are known, then impose an appropriate equation 
of state and also incorporate (jSJ so as to develop a series solution for A and /i. 

It should be noted that solution of the equations (|2l)|) to does not give 
the most general solution of the problem of determining stationary rotating 
fluid sources which match Kerr empty space-time. More general solutions may 
be obtained by appending terms of the type: 

oo 

J2Yk(C)s™ 2k 9 (32) 

k=2 

to each the substitutions (|23|) for U, V, W and H and then to generate a 
further sequence of linear ordinary differential equations for (C) similar to 
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((2T>1) to (|29j). In the absence of an imposed equation of state then the resulting 
equations would be solved subject to the conditions Y& (£ ) = and derivative 
^ (£o) = f° r eacn value of k. The equation of the boundary in case is 
(1231) , When an equation of state is imposed then the equation of the boundary 
is determined from (|24j) and may contain higher harmonic component sin 2fc 9, 
k > 2 depending on the nature of the boundary expressions for {U^} Q and 
{V^} . These expressions will in turn determine the boundary conditions which 
should be imposed on Yj. (£ ) and the derivatives YL (£ ). 

The solution techniques outlined here are of course standard and are there- 
fore not pursued further. The relevance of the analysis to, for example, deter- 
mining potential sources for gravitational radiation or determining more rapidly 
fluid sources matching empty space-time will be examined in future research. 
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